Absorbing states of zero-temperature Glauber dynamics in random networks 
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We study zero-temperature Glauber dynamics for Ising-like spin variable models in quenched 
random networks with random zero-magnetization initial conditions. In particular, we focus on the 
absorbing states of finite systems. While it has quite often been observed that Glauber dynamics lets 
the system be stuck into an absorbing state distinct from its ground state in the thermodynamic 
limit, very little is known about the likelihood of each absorbing state. In order to explore the 
variety of absorbing states, we investigate the probability distribution profile of the active link 
density after saturation as the system size N and (k) vary. As a result, we find that the distribution 
of absorbing states can be split into two self-averaging peaks whose positions are determined by 
(k), one slightly above the ground state and the other farther away. Moreover, we suggest that the 
latter peak accounts for a non-vanishing portion of samples when N goes to infinity while (k) stays 
fixed. Finally, we discuss the possible implications of our results on opinion dynamics models. 

PACS numbers: 05.50.+q, 64.60.De, 75.10.Hk, 89.75.Hc 
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I. INTRODUCTION 

Glauber dynamics [l[ is one of the simplest ways to 
implement the ordering dynamics of Ising spin systems, 
which is defined such that detailed balance holds at equi- 
librium described by the canonical ensemble. When an 
Ising system evades analytical explanation, we can nu- 
merically construct the equilibrium ensemble of the sys- 
tem by running Glauber dynamics until the steady state 
is reached. However, one may ask whether the ensem- 
ble created by Glauber dynamics faithfully represents the 
equilibrium ensemble of the Ising system, a question that 
needs to be addressed case by case with caution. 

The answer is negative when the dynamics occurs at 
zero temperature. Lack of noise and the limited range of 
interaction trap the system in a limited region of phase 
space, or an absorbing state, from which the ground state 
with parallel neighboring spins (the only equilibrium con- 
figuration at zero temperature) cannot be reached. For 
regular lattices, such trappings occur whenever the di- 
mension d is greater than 1. For d — 2, around 30% 
of systems form stripes of domains and do not evolve 
any longer 0, Q. For d > 3, an overwhelming fraction 
of systems fail to reach the ground state. They instead 
stabilize into states of higher energy and complicated do- 
main shapes HQ. The issue of absorbing states has been 
discussed for other kinds of substrates as well, especially 
for complex networks in the context of opinion dynam- 
ics Here, the relevant question is whether a society 



* yj baekOkaist . ac . kr 

t Author to whom all 

msha@chosun.ac.kr 



correspondence should be addressed: 



can reach a consensus through the local majority rule. 
The simplest and most studied case is the Erdos-Renyi 
(ER) random network, which is characterized by the 
system size iV and the average degree (k). Earlier numer- 
ical studies [!, [f| [ll| have shown that ER networks fail to 
reach the ground state in some cases. Combinatorial cal- 
culations by Haggstrom 8] proved a stronger statement 
that the system always fails to reach the ground state in 
the limit when ./V — > oo while (k) is finite. 

Knowing that the system fails to reach the ground 
state, we wonder how close the system can get to it. The 
closeness of the system to the ground state can be repre- 
sented by the active link density I a, which is the fraction 
of links joining spins of opposite signs as illustrated in 
Fig. [1] The ground state has Ia = 0, so greater Ia in- 
dicates that the system is farther away from the ground 
state. While Ia is monotonically decreased by Glauber 
dynamics at zero temperature, it does not evolve any 
more once the system gets stuck into an absorbing state. 
Hence, Ia is a natural choice for characterizing an ab- 
sorbing state. 

Haggstrom || showed that I a in an absorbing state has 
a positive lower bound in the thermodynamic limit, thus 
proving that full relaxation to I a — is impossible. There 
are still more stories to be told about the properties of I a- 
For example, temporal evolution of I a averaged over the 
systems with I a (t) ^ hints at the existence of at least 
two groups of samples with very different behaviors of 
I A (see Fig. [2j where we provide the standard finite-size 
scaling analysis of the quantity measured by Castellano 
et al. [6j]). This leads us to expect that the I a distribution 
has nontrivial structural features worthy of investigation, 
which is still lacking. 

In this paper, we present a systematic analysis of 
the absorbing state distribution of I a produced by zero- 
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FIG. 1: (Color online) An example of a system in which only 
blinkers are flippable. Active links are highlighted in boldface 
(red) and the blinker is marked with a question mark "?". 
Successive flipping of spins at nodes Vi, Vi, . . . , V n leads to 
the merger of two domain walls and further relaxation. How- 
ever, the time required for the process may be long (growing 
like n 2 ), giving a false impression of saturation. This mat- 
ter is addressed by the acceleration algorithm described in 
Sec. EH 



temperature Glauber dynamics in ER random networks. 
In particular, we focus on the influence of parameters 
N and (k) on the distribution. The origin of sample-to- 
sample fluctuations as well as the nature of peaks are also 
discussed. 

The paper is organized as follows. In Sec.|TTl we briefly 
describe zero-temperature Glauber dynamics and define 
quantities of interest. Then, we describe the algorithm 
proposed by Olejarz et al. [H that ensures the saturation 
of the system into absorbing states. Our main results, 
based on extensive numerical simulations, are presented 
in Sec. Mil Finally, we conclude the paper with a sum- 
mary and discussion in Sec. IIVI 



II. 



MODEL 



A. Dynamics and physical quantities 

We consider a system consisting of N nodes. Each 
node is associated with a spin variable cr, (i = 1, . . . , N) 
whose value is either +1 or — 1. Pairs of nodes are 
connected by L randomly distributed links, such that 
(k) = 2L/N is the average degree, or the average num- 
ber of neighbors connected to each node. Let us write 
Aij = 1 if nodes i and j are connected, and = 
otherwise. An active link is a link joining nodes with op- 
positely signed spins as shown in Fig. [T] Then, we can 
define the active link density I a as 



Ia 



(i) 



i<j 



We use zero-magnetization initial states with exactly 
equal numbers of +1 and —1 spins assigned at random. 
At each time step, one spin is randomly chosen and up- 
dated according to zero-temperature Glauber dynamics, 
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FIG. 2: (Color online) Dynamic scaling of (Ia)i a ^o averaged 
over surviving samples that have never reached the ground 
state until time t. Data collapse very well using the scaling 
function (Ia)i A jLo = f(t/ r ), where the relaxation time r scales 
as r ~ In TV. Numerical data are obtained from 10 5 different 
network realizations for (k) — 10 as N varies. The dip struc- 
ture indicates that samples can be divided into two groups: 
the almost fully-ordered case driving the initial fall and the 
quickly absorbed case causing the subsequent rise. 



Eq. That is, a spin flips (changes sign) with the 

following probability, where AIa is the change of Ia re- 
sulting from the flip. 



Flipping probability = 



1, if AIa < 
1/2, if AIa = 
if AIa > 



(2) 



Zero-temperature Glauber dynamics forbids the in- 
crease of I a- Thus, spins aligned with the majority of its 
neighbors are "unflippable." To improve the efficiency of 
our simulations, we employ the rejection-free technique: 
instead of updating one spin among all the N spins, we 
only update one of the flippable spins at each time step. 
The Monte Carlo simulation (MCS) time should then 
be increased by 1/ (number of active spins) per update, 
in order for a fair comparison with the case when the 
rejection-free technique is not used. 



B. Saturation test: absorbing states 

When zero-temperature Glauber dynamics cannot 
lower Ia any further, we say that the system has reached 
an absorbing state. If the system has run out of flip- 
pable spins, it is certainly in an absorbing state since the 
dynamics has stopped. There are, however, absorbing 
states with flippable spins whose neighbors are equally 
split between +1 and — 1. Such spins are called blinkers 
as they repeatedly flip back and forth. The existence of 
blinkers makes the problem of finding an absorbing state 
trickier. 

Blinkers can mislead one to believe that the system has 
saturated to an absorbing state, while I a can still be low- 
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ered as the dynamics continues. Figure [T] illustrates this 
point. It shows a system in which a single blinker is the 
only flippable spin at the moment. The blinker can "hop" 
randomly back and forth between nodes V\ and V n , which 
is effectively a random domain wall motion conserving 
I a- If the random- walk domain wall happens to collide 
with the frozen domain wall on the opposite side of the 
network, the system can go under further relaxation, I a 
eventually going to zero. Hence, the network in Fig. Q] 
has not reached an absorbing state yet. The domain wall 
collision, however, takes some time (proportional to n 2 ) 
due to the randomness of the domain wall motion. Until 
the collision occurs, one may get a false impression that 
I A has saturated and the system has reached an absorb- 
ing state. In order to determine whether an apparent 
saturation of I a truly indicates an absorbing state, we 
employ the algorithm proposed by Olejarz et al. [J] as 
follows: 

1. Until some prescribed MCS time T, let the system 
evolve with zero-temperature Glauber dynamics. 

2. Apply a global infinitesimal field, so that blinkers 
favor one sign over the other and the domain walls 
composed of blinkers are driven in a biased direc- 
tion. The field must not be too strong, lest it should 
make unflippable spins flippable or flippable spins 
unflippable. The field is maintained until domain 
walls collide to decrease I a, or until the system com- 
pletely runs out of flippable spins even in the pres- 
ence of the field. Reverse the global field and do 
the same. 

3. If Ia never decreased in the previous step, the sys- 
tem has saturated to an absorbing state and the 
algorithm terminates. Otherwise, return to the pre- 
vious step. 

To observe the unbiased distribution of absorbing 
states induced purely by zero-temperature Glauber dy- 
namics, one has to determine T such that no lowering of 
I A occurs after applying the acceleration algorithm. This 
is quite a time-consuming task. So we just accept the 
absorbing state found by the acceleration algorithm even 
if the algorithm has lowered I a- This inevitably biases 
the absorbing state distribution of I a, but we have made 
sure that the change in the distribution is negligible if T 
is larger enough than r. Here, we set T = 100 In N. 

III. NUMERICAL RESULTS 

Utilizing the acceleration algorithm to find true ab- 
sorbing states of zero-temperature Glauber dynamics in 
ER networks, we obtain the probability distribution func- 
tion P(Ia) of the active link density I a in absorbing 
states through extensive numerical simulations covering 
both different runs on a given ER network and different 
network realizations. P(Ia) is typically bimodal with 
the left peak near Ia = and the right peak at a larger 
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FIG. 3: (Color online) Probability distribution functions of I a 
in absorbing states (a) at (k) = 1, 2, 3, 4 with N — 2 10 and 
(b) at TV = 2 9 , 2 11 , 2 13 with (k) = 3. Numerical data are ob- 
tained from 10 2 zero-magnetization initial configurations and 
10 network realizations. The Ia distribution is essentially bi- 
modal, which implies that it is a mixture of two components 
that need to be analyzed separately. 

value of I a (see Fig. [3]). Such a shape of the distribution 
is consistent with the temporal behavior of (Ia)i a ^o m 
Fig. [2j both imply that there are two groups of samples 
with very different properties. From this perspective, 
we analyze P (I a) assuming that the distribution always 
consists of two components, the left and the right. P (Ia) 
may appear unimodal in some cases, as illustrated by the 
curves at (k) = 1 and (k) = 2 in Fig. |3fa). Even in such 
cases, we can assume that there is an unobserved peak 
of zero (or extremely small) magnitude. An example is 
shown in Fig. [3Jb) , where the left peak seems to shrink 
indefinitely as N grows, leaving only the right peak to be 
observed. 

In the following subsections, we systematically inves- 
tigate the finite-size effect on P (I a) with the goal of ob- 
taining the shape of the distribution in the thermody- 
namic limit. We present numerical evidence that the Ia 
distribution consists of at most two delta peaks in the 
thermodynamic limit, and that a nonvanishing propor- 
tion of samples belongs to the right peak. Clues as to the 
nature of the peaks and the origin of sample-to-sample 
fluctuations of Ia are also discussed. Finally, we provide 
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FIG. 4: (Color online) The fraction of samples belonging to 
the right component, 4>r, is calculated from the I a distribu- 
tion, (a) If TV is fixed, the value of (j>R decreases with (k). (b) 
If (k) is fixed, <f)R eventually increases with TV. Larger TV is 
required at greater (fc) to observe the increasing trend. The 
lines are a guide to the eyes. 

evidence that P {I a) may have only a single delta peak 
far from I a = in the thermodynamic limit. 

A. Proportions of components 

We begin our analysis with the finite-size effect on the 
fraction 4>r of samples belonging to the right component 
of the I a distribution. The component corresponds to 
samples stuck far away from the ground state, so higher 
indicates greater difficulties in reaching the ground 
state. We are also interested in whether the I a distribu- 
tion is still bimodal in the asymptotic limit. 

Measurement of 4>r requires the criterion for distinc- 
tion between the two components of the I a distribution. 
Without loss of generality, we choose the local minimum 
point Z m i n between the maximum point of the left peak 
^max and the maximum point of the right peak Z^ ax as the 
border between two peaks. Therefore, 4>r is the fraction 
of samples ending up on the right side of Z m i n . If there 
are multiple candidates for Z m i n , the largest (rightmost) 
one is chosen for a conservative estimation of (f>ji . 

According to Fig. 0|a), (f>R decreases monotonically 



with (k) when TV is fixed. This is a natural consequence 
of the fact that nodes are better informed of the global 
situation when the connectivity is greater. On the other 
hand, 4>r has a nonmonotonic dependence on TV when (k) 
is fixed, as in Fig.[4jb). <pR eventually increases with TV, 
which implies that the right-peak component accounts 
for a nonvanishing portion of samples in the thermody- 
namic limit. We note that the crossover occurs at greater 
TV as (fc) gets larger, so observing highly connected net- 
works of small sizes may mislead one to conclude that 
the right component becomes negligible in the thermody- 
namic limit. The behavior of 4>r at greater TV shows that 
this is not likely to be the case. On the contrary, one may 
even claim from our result that the right peak is actu- 
ally the only observed peak in the thermodynamic limit, 
as can be numerically confirmed for sufficiently small (k). 
We shall get back to this bold claim later, providing more 
numerical evidences. For now, we simply interpret the 
numerical observation for <f>n as evidence that the right 
component is not an artifact due to the finite-size effect. 



B. Means and widths of peaks 

Now we take a closer look at the shape of P {I a), and 
measure the finite-size effect on positions and widths of 
the two peaks, which are defined as follows. Three ex- 
trema (l^^, ^Sax; ^min) defined in the previous subsec- 
tion divide the domain of I a into four intervals. The 
medians of the intervals are denoted by l\ < h < h < h 
in the ascending order. Then, the width of the left (right) 
peak is defined as wl = h~h {wr = Z4 — Z3), and the po- 
sition is represented by the mean value [(^a)]/, ([(Ia)]r) 
which is the average of I a among samples belonging to 
the interval between l\ and I2 (h and The behaviors 
of those peak properties are plotted in Fig. [5j We note 
that both peaks are self-averaging: the mean values show 
signs of saturation, and the widths seem to decay asymp- 
totically like TV -1 / 2 . Hence, in the thermodynamic limit, 
we claim that almost all the samples belong to two sharp 
peaks whose positions are determined solely by (k). Here 
follows detailed descriptions for the analysis of each peak. 

1. Analysis of the left peak 

We present the left-peak properties for a rather limited 
range 3.25 < (k) < 4.25, due to the following limitations 
of our study. The left peak may not be observed for 
small values of (k) as shown in Figs. 02a) and |3J so its 
properties cannot be measured in such cases. On the 
other hand, the left peak is too sharp and too close to 
zero for large values of (k) , so we hit the resolution limit 
of peak position and width. In such cases, the left-peak 
properties show the trivial behavior of a peak of unit 
width (proportional to 1/N(k}) located at I a = 0. 

Figure IHJa) shows that the value of decreases 
with (fc), and approaches positive limits as TV grows. This 
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FIG. 5: (Color online) The means and widths of two peaks of P{Ia) are measured, (a) The left-peak mean [{Ia)] l is plotted at 
(k) = 3.25, . . . , 4.25 increased by 0.25 from top to bottom, (b) The left-peak width w L is plotted for the same parameter sets as 
(a), (c) The right-peak mean [{/a}]^ is plotted for (k) = 1, . . . , 7 increased by 1 from bottom to top. (d) The right-peak width 
w a is plotted for the same parameter sets as (c). Note that the curves at (fc) = 1, 2 have been vertically shifted by respective 
factors of 4 and 5 to avoid overlapping with other data. All the data are obtained from 10 3 different initial configurations and 
10 3 network realizations. The lines are guide to the eyes. 



is consistent with the lower bound of I a in the thermo- 
dynamic limit derived by Haggstrom [8| which is positive 
and decreases with (fc) in the parameter range observed 
in our simulations. The vertical drops of I a indicate that 
the peak becomes too sharp and too close to zero at 
smaller values of N, so becomes indistinguishable 

from zero in those ranges. 



Figure Ob) shows that the value of wt, has a crossover 
between two scaling regimes: N- 1 and TV" 1 / 2 . As 
pointed out in previous paragraphs, the left peak is so 
narrow for small values of A~ that its width is just given 
by the unit of I a, which is inversely proportional to A. 
As A~ gets larger, the unit of Ia decays faster than the 
true peak width, revealing the A^ -1 / 2 dependence. This 
crossover to the true scaling regime occurs for the larger 
values of (k) as A gets larger, which is obviously related 
to the shifting of vertical lines in Fig.[5ja), implying that 
the two artifacts are of the same origin. 



2. Analysis of the right peak 

The right-peak properties are presented for a wider 
range of parameters, as the limitations present in the case 
of left-peak properties are not as severe. For the sake of 
clarity, we have reduced the resolution of (k) compared 
to the case of the left peak. 

While Fig. [5jc) clearly shows that [(Ia)]r saturates as 
A~ grows and increases with (k), wr exhibits more com- 
plicated behaviors as shown in Fig. [5fd). For (k) > 3, 
the curves seem to follow the N^ 1 / 2 scaling at small N, 
but deviate from it as A" grows. On the other hand, for 
(k) < 3, the curves are increasingly well-fitted by A^ 1 / 2 
as A" increases. We claim that wr decays somewhat like 
A -1 / 2 , but deviates from it around the regime where the 
left peak is just about to disappear (or, equivalently, just 
about to appear). The data for (k) = 3 in Fig. [5fd) 
and the distribution shapes in Fig. [3fb) give a nice illus- 
tration of our claim. When A is small, the two peaks 
are clearly distinguishable, so the A -1 / 2 behavior is eas- 
ily observed. As N grows, the left peak dwindles and 
the left component becomes indistinguishable from the 
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right component of the I a distribution. This makes wr 
appear broader than it really is as we include part of 
the left component in the right component by mistake, 
hence wr appears to decay slower than TV -1 / 2 . As N 
grows more, the left peak becomes completely unobserv- 
able, causing the rapid width decay and the subsequent 
return to the TV" 1 / 2 behavior. If we had observed the 
system at greater N, the curve for (k) = 7 would also 
have shown the rapid width decay and the subsequent 
return to the iV -1 / 2 behavior as was shown by the curve 
for (k) = 3. 

Taking all those considerations into account, we con- 
clude that both components of P (Ia) are self-averaging. 



C. Nature of the peaks 

Based on the results presented so far, we claim that 
the I a distribution in the thermodynamic limit can be 
represented as a combination of at most two delta peaks. 
The nature of the left peak is easily understood: it rep- 
resents the samples that almost reach the ground state, 
with a small (but still finite) fraction of active links con- 
necting dangling parts of the system. However, we pose 
the following question: why should there be another well- 
defined peak at a much larger value of Za? 

In order to check whether the ER network is the only 
substrate with the observed right-peak behavior, we also 
test zero-temperature Glauber dynamics in the regular 
random network where all nodes have the same degree k. 
The results are shown in Fig. [6ja) . It is remarkable that 
zero-temperature Glauber dynamics in even-degree reg- 
ular random networks never produces any "right peak" 
(only a single peak is formed at or around I a = 0, which 
may be called the "left peak"), which is why the data 
are not plotted at even values of (k). This seems to be 
due to large fluctuations provided by blinkers that can 
appear only in even-degree nodes. On the other hand, 
odd-degree regular random networks produce results very 
similar to ER networks. Lack of degree fluctuations in 
odd-degree regular random networks yields only small 
vertical offsets of [(Ia)]r, not significantly altering the 
right-peak properties observed in ER networks. 

To get some hints as to whether the location of [(Ia)]r 
has any physical meaning, we compare [{Ia)]r obtained 
from the previous simulations in ER networks (whose 
initial magnetizations are zero, mo = 0) with the same 
quantity obtained from initial states characterized by the 
initial active link density Iq. For the sake of convenience, 
let us denote the latter by a new notation [(Za)]^- Since 
I A always decreases with time in Glauber dynamics, sam- 
ples initialized by Iq cannot produce right peaks at values 
of Ia larger than Iq. Therefore, if [(Ia)]r observed in the 
previous results are truly special, we expect a drastic dif- 
ference in the behavior of {(Ia)] 1 r between the samples 
with Z > [(Ia)]r and those with Iq < [(Ia)]r since the 
former can produce a peak around [(Za)]r while the latter 
cannot. 



Before moving on to the discussion of results, a re- 
mark on the preparation of samples with Iq is in order. 
Those samples are constructed from zero-magnetization 
spin configurations by random flippings of spins that de- 
crease \Ia — Iq\- It should be noted that this method does 
not scan the entire configuration space of the samples 
characterized by Zo, since the magnetization of each sam- 
ple prepared by the method is likely to be very close to 
zero. However, since the samples belonging to the right 
component have magnetization close to zero in the ab- 
sorbing states, our initialization method still serves the 
purpose of generating samples that are similar to the 
right-component samples but different only in the value 
of l A - 

Figure [D(b) indeed confirms our previously described 
expectation: the samples with Z > [(Ia)]r have [(Za}]# = 
[(Za)].r, while the right peak disappears in the samples 
with Iq < [{Ia)]r (the right peak may still remain for 
a few more steps of (k) as the right peak produced by 
samples with too = does have a tail to the left of [(Za)]#, 
which is smaller than Iq). This indicates that the value 
of [(Za)^ is indeed physically meaningful since it is the 
only value of active link density that can be sustained 
when the magnetization of the system is close to zero. 



D. Origin of sample-to-sample fluctuations 

Now we discuss how the fate of each sample is deter- 
mined, or what determines whether the system evolves 
into the left or the right component of the I a distribution. 
Such sample-to-sample fluctuations come from three fac- 
tors: different network structures, different initial config- 
urations, and the stochasticity of dynamics. 

We separate the effect of network structure from that 
of the other factors. A natural parameter describing the 
broadness of the Ia would be its variance, which we may 
define in the following two ways: 

° 2 = [<&>] -Wa)} 2 , 
°' 2 = [<&> - (Ia) 2 } ■ 

Here, [. . .] represents an average over different network re- 
alizations, and (. . .) represents an average over the other 
factors sharing the common network structure. There- 
fore, a 2 is the variance of Ia that does not distinguish 
between contributions from the three factors mentioned 
above. On the other hand, a' 2 is the different variance of 
Ia due to different initial configuration and the stochas- 
ticity of dynamics, averaged over different network real- 
izations in the final step of calculation. The difference 
between a 2 and a' 2 can be interpreted as the fluctua- 
tions of I a solely due to differences in the network struc- 
ture. Strictly speaking, a 2 = a' 2 implies that (I a) is the 
same, irrespective of the exact structure of ER networks 
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FIG. 6: (Color online) (a) The right-peak mean of Ia measured for 10 7 ER network samples (lower points) and 10 6 regular 
random network samples (upper points). Lack of degree fluctuations in odd-degree regular random networks does not make 
much difference in the right-peak mean, while the right peak completely disappears in even-degree regular random networks, 
(b) The same quantity measured for different initial magnetization mo or initial active link density Iq in fO 6 ER network 
samples. The result suggests that for each value of (k), the right peak corresponds to the only stable spin configurations with 
approximately zero magnetization. 
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FIG. 7: (Color online) Samplet-to-sample fluctuations caused by the contributions of network structure and initial spin con- 
figuration are measured for 10 2 initializations of f 5 network realizations, (a) Variances of Ia measured in two different ways 
(refer to the main text for details) agree well with each other, implying that the mean of Ia distribution is unaffected by the 
exact structure of each particular ER network, (b) A similar result is observed for Binder cumulants of I a- Here open symbols 
represent (er 2 , Ui) and solid lines represent (o 72 , U4), respectively. 



as proven below, 



moments of I a, such as Binder cumulants Ua and U. 
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a' 2 -a 2 = [<^> 2 ] - [(l A )Y 
= '((Ia) [(l A )]f 
(Ia) = [(U)] ■ 



Figure [Zta) shows that a 2 and a' 2 are close to each 
other for various combinations of N and (k). Therefore, 
the network structure does not contribute much to the 
mean values of I a- The same can be said for the higher 
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which are shown in Fig. [7f b) . Hence, we can conclude 
that the shape of the I a distribution is not really changed 
by differences in the network structure. Similar results 
are observed when we distinguish between the contribu- 
tions of the initial configuration and the stochasticity of 
dynamics. Again, we can say that the shape of the I a dis- 
tribution is not affected by the difference among initial 
configurations. 
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FIG. 8: (Color online) The ./V-dependence of three specific 
values of (k): ku (o) , k a (°) , ki (x). They are possible 
indicators of the disappearance of the left peak (see the main 
text for their precise definitions). 



E. P (I a) in the thermodynamic limit 

Finally, we return to the question posed early in this 
section: is the right peak the only lasting peak in the 
thermodynamic limit? To answer this question, we de- 
fine possible indicators of the point at which the left peak 
becomes unobservable. The TV-dependence of those indi- 
cators may reveal whether the left peak survives in the 
thermodynamic limit. 

The curves shown in Fig. [7] reflect the (A;)- dependence 
of P {I a)- As (k) grows, we observe the shift of the single 
(right) peak to the right, and then the formation of the 
left peak, and finally the increasing inter-peak separation 
with the dominance of the left peak, as Fig. E^a) illus- 
trates. This produces the single-peaked curves as shown 
in Fig. [7ja). Although not exhibited here, [(Ia)] shows 
a similar behavior. Meanwhile, the Binder cumulant U4 
indicates the separation between two peaks, whose value 
decreases far below zero as the separation increases in- 
definitely with (fc). 

To investigate the finite-size scaling of Fig. we define 
the three specific values of (k) as 

• kjj: the largest (k) that satisfies U4 = 0, 

• k a : the location of max a 2 , 

• ki: the location of max [{Ia)]- 

Taking those three points as possible indicators of the 
critical value of (k) at which the left peak is barely ob- 
servable, we observe that all the three values increases 
indefinitely with N (see Fig. |8]). This implies that as N 
goes to infinity, the I a distribution is single-peaked at any 
finite (k), and that the peak is far away (as it is actually 
the "right" peak) from l A = 0. Hence, this observation 
again supports our previous claim that the right peak 
is actually the only lasting peak in the thermodynamic 
limit. 



IV. SUMMARY AND DISCUSSIONS 

In summary, we have studied the absorbing states of 
zero-temperature Glauber dynamics in quenched ER- 
type random networks with zero-magnetization initial 
states. Each absorbing state has been represented by its 
active link density I a- We have found that the structure 
of the Ia distribution consists of two self-averaging peaks, 
the left one with samples almost reaching the ground 
state and the right one with samples divided by large 
stable domain walls. The mean of the left-peak (right- 
peak) component seems to be a monotonically decreasing 
(increasing) function of (k). In particular, we focus on 
the finite-size scaling analysis of various physical prop- 
erties for P(Ia), which suggests that the fraction of the 
right-peak component stays finite in the thermodynamic 
limit, even if the value of (k) is pretty high. We have 
also gathered circumstantial evidence that the right peak 
is actually the only lasting peak in the thermodynamic 
limit. 

Our results also have some implications on opinion dy- 
namics studies. Some earlier studies of zero-temperature 
Glauber dynamics as an opinion dynamics model have 
simply assumed that when there are sufficiently many 
links, the portion of samples that do not reach the ground 
state becomes negligible. Thus, they measured "consen- 
sus time," which is the time required for the sample to 
reach the ground state in order to obtain the scaling rela- 
tionship between the system size and the speed of relax- 
ation. We suggest that consensus time is not a suitable 
measure since almost all the samples fail to reach the 
ground state in the thermodynamic limit and increas- 
ingly many of them end up far away from the ground 
state as (fc) gets larger [lSj. Data collapse shown in 
Fig. [2] (the results of Q as well) is a better numerical ap- 
proach to such a problem. The final location of the right 
peak is very well-defined and it seems to represent the 
only value of Ia that yields stable up-down symmetric 
configurations. Therefore, interfacial noise that causes 
sufficiently large fluctuations in Ia always drives the sys- 
tem out of absorbing states and helps them relax to the 
ground state, as illustrated by the role played by blinkers 
in even-degree regular random networks. This leads us 
to expect that there should be an analytical argument 
for the location of the right peak, which we unfortu- 
nately could not formulate. As the right peak observed 
in odd-degree regular random networks shows a quali- 
tatively similar behavior, an intuitive argument for the 
nature of the right peak in one substrate may help us 
understand the corresponding phenomenon in the other 
substrate. 

Finally, we have shown that the properties of P(Ia) 
are not much affected by either the particular network 
structure or initial configuration, where both peaks can 
be reached purely from the stochasticity of dynamics. 
This may be due to the characteristic homogeneity of 
ER networks, as any randomly generated ER network 
in any randomly chosen symmetric initial configurations 
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is likely to be very similar to another ER network pre- 
pared in the same way. It would be interesting to check 
whether the same observation can be made for more topo- 
logically heterogeneous substrates such as scale-free net- 
works. One should note that this argument for the "irrel- 
evance of structure and configuration" should be taken 
carefully. As Uchida and Shirayama have shown the 
biased choice of initial configuration does affect the prop- 
erties of P(Ia)- The contributions from such exceptional 
cases are negligible in the unbiased ensemble of samples 
we have considered in this paper. 
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